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Abstract 

The isovector dipole El strength in 92,94,96,98, loOjy^Q jg analyzed within the self-consistent separable 
random-phase approximation (SRPA) model with Skyrme forces SkT6, SkM*, SLy6, and SkI3. The special 
attention is paid to the low-energy region near the particle thresholds (4-12 McV), which is important for 
understanding of astrophysical processes. We show that, due to a compensation effect, the influence of 
nuclear deformation on El strength below 10-12 MeV is quite modest. At the same time, in agreement with 
previous predictions, the deformation increases the strength at higher energy. At 4-8 MeV the strength is 
mainly determined by the tail of El giant resonance. The four Skyrme forces differ in description of the 
whole giant resonance but give rather similar results below 12 MeV. 



* Permanent address: Tech. Universitat Dresden, Inst, fiir Analysis, D-01062, Dresden, Germany, 
t Present address: Department of Physics, University of Notre Dame, Indiana, 46556, USA. 



I. INTRODUCTION 



The isovector giant dipole resonance (GDR) remains to be a subject of intense study, now with 
the accent to exotic nuclei and astrophysical problems m. Besides the GDR is actively used for 
inspection and upgrade of the modern self-consistent mean-field approaches [2jl3|], e.g. those based 
on Skyrme forces [J, . 

In this paper we will investigate GDR in a chain of even-mass isotopes ^^"^'''^Mo. This chain 
is of particular interest because here we have photoabsorption experimental data not only above 
but also below the particle emission thresholds, down to 4 MeV [a, Ql- Being rather weak, the El 
strength near the thresholds is however important for understanding astrophysical processes, e.g. 
of the stellar photodisintergation rate 0,0, [3]. In principle, this strength should depend on such 
factors as nuclear size and deformation. Since the chain ^■^"^'''^Mo involves spherical (A=92,94,96) 
and deformed (A=98,100) nuclei 1]|, it is suitable for estimation of both factors. 

The El strength in the molybdenum chain has been recently explored in the random-phase- 
approximation (RPA) with the phenomenological Nilsson [12] and Woods-Saxon j?] single-particle 
potentials. In both cases, the important role of the deformation was found. Namely, it was 
predicted that prolate/triaxial deformation significantly increases the El strength at 10-14 MeV. 
However, the calculations are not self-consistent and so an additional analysis within more 

involved microscopic models is desirable. In the present paper we explore the dipole strength 
in ^^~^°^Mo in the framework of the self-consistent separable RPA (SRPA) with Skyrme forces 
SRPA covers both spherical 1^ and deformed 15, Q, [lit and so is the proper tool for 



the present analysis. Factorization of the residual interaction minimizes the computational effort 
and therefore allows the systematic exploration. The model was already used for the analysis of the 
GDR in different mass regions, including drip line and superheavy nuclei [l7 . 3]. Note that unlike 



the previous Skyrme- RPA studies of the low-energy El strength [10|, where the GDR deformation 
splitting was introduced phenomenologically, SRPA treatment of the deformation effects is fully 
self-consistent. 

We will show that conclusions on the impact of nuclear deformation should be amended 

in the sense that the result strictly depends on the energy region. Namely, the deformation indeed 
increases the El strength at E > 10-12 MeV but, at the same time, has a minor impact at lower 
energy, i.e. near the particle thresholds. The later is caused by a compensation effect of the GDR 
branches. This effect is quite general and becomes apparent just in the low-energy regions of 
astrophysical interest. 
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II. CALCULATION SCHEME AND GROUND STATE PROPERTIES 



The calculations are performed within SRPA 

n " ~ 

forces, SkT6 [19], SkM* [20], SLy6 



with the representative set of Skyrme 



22| . These forces expose a variety of features 
17[. Amongst them, the force SLy6 provides 



21[, and SkI3 
(effective masses, etc) relevant for the GDR 

m 

the beest compromise for the description of the GDR in heavy nuclei 
residual interaction involves contributions from the time-even densities (nucleon p, kinetic r and 
spin-orbital 9), time-odd current j, direct and exchange Coulomb terms, and pairing (with delta 



15|, Q, Q. The SRPA 



forces at the BCS level) 



13|, 
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FIG. 1: Energy surfaces (= binding energies) in Mo isotopes, obtained in quadrupole-constrained mean-field 
calculations with Skyrme forces SkT6, SkM*, SLy6, and SkI3. The equilibrium axial deformations (3 are 
indicated by arrows. 
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TABLE I: Experimental thresholds for (7,n), (7,p), {j,2n), (j,np), and (7, 2p) reactions 23[ and equilib- 
rium deformations P calculated with the forces SkT6, SkM*, SLy6, and Skl3 in ^'^~^°°Mo. 



A 


Thresholds (MeV) 


f3 


En 


Ep 


E2n 




E2p 


SkT6 


SLy6 


SkM* 


Skl3 


92 


12.7 


7.5 


22.8 


19.5 


12.6 


0.0 


0.0 


0.0 


0.0 


94 


9.7 


8.5 


17.7 


17.3 


14.5 


0.0 


0.0 


0.0 


0.0 


96 


9.2 


9.3 


16.5 


17.8 


16.1 


0.0 


0.05 


0.0 


0.18 


98 


8.6 


9.8 


15.5 


17.9 


17.2 


0.0 


0.22 


0.0 


0.24 


100 


8.3 


10.1 


14.2 


18.0 


19.5 


0.0 


0.27 


0.0 


0.24 



The axial equilibrium quadrupole deformation 



/3 



1 



5 Z <r^ >r 



drpp{r)r'^Y2o 



(1) 



(where Pp{r) is the proton density in the ground state and < >p= J drpp{r)r'^ / Z is the r.m.s. 
proton radius) is determined by minimization of the total energy, see Fig. 1 and Table HI As 
is seen from Fig. 1, the nuclei ^^"^'''^Mo are soft to j3, especially the heavy isotopes. Partly 
this is because of their transition character. Besides, ^^Mo and ^"^Mo are probably triaxial (with 
e2 = 0.18, 7 = 37° and €2 = 0.21, 7 = 32°, respectively {ll|, [3]). In the present study the triaxiality 
is omitted, which may also lead to a wide plateau and local minima for /3 < 0. Altogether, there 
is appreciable ambiguity in determination of the equilibrium deformation j3. In ^^'^'^'^Mo we see a 
significant variation of (3 with the Skyrme force. While SkT6 and SkM* favor a spherical shape 
(in contradicts with significant 62 in Refs. Ul, [l^])) the forces SLy6 and SkI3 give more reasonable 



results. So just SLy6 and SkI3 will be mainly used in the further analysis of the deformation 
effects. Note, that triaxiality can cause an additionally spread of the El strength. In the present 
study this effect is masked by the proper Lorentz smoothing of the strength. 



The photoabsorption (in mb) is computed as [24 1 



where 



o^{E) = 4.01 E S{E) 



At=0,l V 



(2) 



(3) 



is the strength function with the Lorentz weig ht C = A/(27r[(£; - E^f + A2/4]) and isovector 



transition operator fEi^, = I]p=i ^P^i/i(^p) " En=i • ®' ^'^'^ 

mark eigenvalue and eigenfunction of z^-th RPA state, and ^'o is the ground state eigenfunction. 
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The summation runs over the RPA spectrum until Ecut = 45 MeV. The Lorentz function with 
the averaging parameter A is used to simulate the broadening effects beyond SRPA (triaxiality, 
escape widths, and coupling with complex configurations). Besides, this smoothing allows to avoid 
unnecessary details of the calculated strength which in any case are not resolved in experiment. 
Following 



15 . [itI . [isl . the averaging A=2 MeV is optimal. SRPA allows a direct computation 
of the strength function (without finding the manifold of RPA states), which greatly reduces the 
effort. 

The calculations use a large basis space of single-particle states, from the bottom of the potential 
well up to ~ +16 MeV. The integral photoabsorption 

"Ecut 

dE a^{E) (4) 







exhausts up to 98% of the energy-weighted sum rule EWSR = 9(/i.e) /(Svrm^^) • NZ/A wit 



1 the 



17 



251]. 



isovector effective mass m* arising due to the velocity-dependent densities r, 9, and j 
So our basis is indeed large enough to explore the GDR. The EWSR should not be confused with 
the Thomas-Reiche-Kuhn (TRK) sum rule Titrk = 60NZ/A mb MeV which deals with the bare 
nucleon mass m. 

III. RESULTS AND DISCUSSION 

Results of the calculations for El strength are given in Figs. 2-7. 

In Figs. 2 the photoabsorption in semimagic spherical ^^Mo is presented. The logarithmic scale 
is used for a convenient comparison with the experimental data. Different smoothing parameters 
A are used. It is seen that the most appropriate agreement with experiment takes place for A=l 
and 2 MeV, especially in the low-energy 4-12 MeV of our interest (the same for other isotopes 



94-i00]yfQ^_ Taking also into account our previous results for the GDR 15, [itI . [isl . we chose for 
the further analysis A=2 MeV. As is seen from Fig. 2, for this averaging we generally reproduce 
the GDR energy and width. At the same time, the calculated cr^{E) systematically exceeds the 
experimental data at the GDR top and falls short at its right flank. This discrepancies can be 
partly caused by omitting the complex configurations. Being included, those configurations could 
redistribute the strength from the middle to the right side of the GDR, thus improving agreement 
with the experiment. Besides, the discrepancy at the GDR top can be also an artifact of the 
comparison of the calculated photoabsorbtion with the separate (7,^*) and (7, xn) channels. 

In the logarithmic scale the photoabsorption obtained for different Skyrme forces looks rather 
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FIG. 2: Photoabsorption (in logarithmic scale) in ^^Mo calculated with Skyrme forces SkT6, SkM*, SLy6, 
and SkI3 with the averaging A=0.25 MeV (dotted curve), 1 MeV (solid curve) and 2 MeV (bold solid curve). 
The experimental (7, 7'), (7,^), and (7, xn) data j^, 7| are given by boxes, triangles, and stars in the energy 
intervals 5.1-12.5, 7.9-20.3, and 12.4-29.7 MeV, respectively. 
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FIG. 3: The same as in Fig. 2 but for linear scale and averaging A— 2 MeV only. The dash curve exhibits 
the unperturbed two-quasiparticlc (2qp) results. 



similar. To distinguish the difference between predictions of different forces, one should switch to 
the linear scale, which is done in Figs. 3 and 4 for semimagic ^^Mo and deformed ^'^'^Mo. Now 
it is seen that in spherical ^^Mo the best description is for SLy6. For deformed ^'''^Mo only SLy6 
and SkI3 results are presented. The forces SkT6 and SkM* predict for ^^'^Mo the spherical shape 
and so obviously fail here. The calculations for ^'^'^Mo overestimate the experimental data (as 
photoabsorption over (7,xn)) but keep EWSR and provide an acceptable agreement for the GDR 
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FIG. 4: The same as in Fig. 3 for ^"Mo and forces SLy6 and SkI3. The experimental data 



for (7, 7') 



at 4.1-8.1 MeV and {■^,xn) at 8.3-27 MeV are given by boxes and stars, respectively. The GDR branches 
fi — (small bump) and /i = 1 (twice larger bump) exhibited by the dash curve demonstrate the deformation 
splitting of the resonance. 



energy centroid and width. 

As was mentioned in the introduction, the El strength near the particle thresholds is of par- 
ticular interest for some 

understand the main physical mechanisms responsible for the evolution of this strength with the 
neutron number N in the isotope chain. First of all, this evolution is determined by the empirical 
rule Egdr = SlA'""/^ MeV [2^ relating the GDR energy and nuclear mass number A. The higher 
N (and so A), the more the GDR downshift and stronger the GDR tail (and relevant El strength) 
near the thresholds. This size factor seems to dominate. However, two other effects, nuclear de- 
formation and internal El strength in the region (pygmy resonance), can also come to play and 
influence the general size trend. 

These effects are inspected in Figs. 5-7. The impact of deformation is illustrated in Fig. 5, where 
the El strength calculated at zero and non-zero deformations is compared with the experimental 
data y, Q]- We get the nonzero deformation only in ^^Mo and ^'^'^Mo for the forces SLy6 and 
SkI3. In all other cases the calculations give /3 =0, see Table 1. Fig. 5 shows that for /3=0 all the 
forces provide more or less acceptable agreement with the experiment, the best for SkT6 and SkI3. 
Inclusion of the deformation changes the results. While at E' > 12 MeV we have, in accordance 
with Refs. [3, increasing El strength, in the interval of our main interest, E < 12 MeV, we 
see its modest decrease. The effect is prominent and even stronger than the change in El strength 
for the neighbor isotopes. 

The deformation effect for E < 12 MeV can be explained (for both prolate and oblate shapes) 
by destructive competition of two factors, deformation shifts of = and = 1 branches of the 
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FIG. 5: The calculated a)-d) and experimental ^ 3] e) integral low-energy photoabsorption in Mo isotopes. 
The calculated results are exhibited for non-zero (= equilibrium in ^^'^^^Mo for SLy6 and SkI3) and zero 
{— equilibrium in other cases) deformations by solid and dotted curves, respectively. In all the panels the 
sequence of curves for 92,94,96,98, loojy^Q j^^^g ^j^g same order: from the lowest for ^^Mo to the highest for 
loOMo. 



GDR. Indeed, the deformation shifts the GDR branches in opposite directions thus minimizing 
the resulting deformation impact. Though the branch /j, =1 is twice stronger than /u =0 one, its 
deformation shift is less, hence a strong mutual compensation of both ^ = and = 1 deformation 
impacts. Being strong, the compensation is usually not complete. So, depending on the concrete 
case, one may finally observe a modest decrease or increase of the low-energy El strength. Of 
course, the compensation holds at the energies far enough from the lower GDR branch (which is 
just the case for E < 12 MeV). Instead, while approaching the lower GDR branch, we always gain 
El strength with the deformation, as was earlier found in Refs. 

These arguments are illustrated in Fig. 6 for prolate ^'^'^Mo. The calculations were done for SLy6 
with averaging parameters A = 2 and 1 MeV. It is seen that for A = 2 the deformation-induced 
/i = and 1 contributions indeed strictly compensate each other, thus leading to a slight decrease 
of El strength at 9 < < 12 MeV and no effect at the lower energy. Instead, in accordance with 
, the strength grows at ii^ > 12 MeV. The minor deformation effect at 9 < ii^ < 12 MeV 
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E [MeV] E [MeV] 

FIG. 6: The low-energy photoabsorption in ^""Mo, calculated with the force SLy6 for f3 —0 (dotted curves) 
and equilibrium deformation [3 —0.27 (solid curves). The averaging A =2 MeV (left) and 1 MeV (right) is 
used. The total strength as well as the strengths of the branches /x = and = 1 are indicated by arrows. 
The experimental data Q| are depicted by full boxes. 




E [MeV] E [MeV] 

FIG. 7: The low-energy photoabsorption in deformed ^^Mo, calculated for the forces SLy6 and SkI3 with 
(solid curve) and without (dash curve) the residual interaction. To see the fine structure, the small averaging 
A =0.1 MeV is used. The unperturbed photoabsorption is decreased 10 times as indicated. The vertical 
arrows show (7,n) and (7,p) thresholds. 

becomes invisible for the less averaging A = 1 MeV, which is explained by weakening the GDR 
low-energy tail. Note that for A = 2 MeV we get much better agreement with the experiment than 
for A = 1 MeV, which additionally justifies the large averaging as the best choice. 

Finally, Fig. 7 demonstrates for deformed ^^Mo the role of internal El excitations in forming 
the low-energy strength. To make visible a fine structure, the small averaging A =0.1 MeV is 
used. It is seen that for both forces the El strength between 4 and 7-8 MeV does not exhibit any 
internal structure while the structure at 7-9 MeV is weak. So, the El strength for E < 9 MeV 
is mainly determined by the GDR tail (the same for other Mo isotopes). This means that all the 
effects discussed above for the GDR are indeed relevant for the energies near and below the particle 
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thresholds. The absence of El structures at 4-8 MeV is natural since 2qp excitations with AA^=1 
{Af is a principle shell quantum number) lie at a higher energy, see Figs. 3 and 4. Note that for 
£ > 9- 10 MeV the El^*eng.h H.e. and exhibits more appreciable st™cU.e (perhaps the py.my 
resonance observed in jQ])- 



IV. SUMMARY 



The El strength in 92,94,96,98,i00]y[Q jg investigated in the framework of the self-consistent sepa- 
rable RPA method for the set of Skyrme forces SkT6, SkM*, SLy6, and SkI3. The main attention 
is paid to low-energy strength below the particle thresholds, which is of a keen interest for astro- 
physical problems [a, Q]. To our knowledge, this is the first Skyrme- RPA study of El strength 
in Mo isotopes with the self-consistent treatment of deformation effects. Some important factors 
(triaxiality, coupling with complex configurations, escape widths) are simulated by the Lorentz 
smoothing and others (pairing impact, energy dependence of the smoothing) need an additional 
analysis, which makes our description of El strength yet tentative. Nevertheless, some useful 
conclusions can be done. 

Q. 



We confirmed our previous statement 



15 



13, 



la ] that the force SLy6 with the Lorentz 



smoothing A = 2 MeV gives the most reasonable description of the whole GDR. The low-energy 
El strength is shown to be mainly determined by the GDR tail. In spherical 92,94,96]y[Q^ ^^lis 
strength is well described by all the forces. 

It is found that the deformation impact in the low-energy El strength depends on the particular 
energy interval. While approaching the GDR, E > 12 MeV, we get a definite increment of the 
strength with the deformation (as was earlier found in Refs. ^jll^])- However, at -B < 12 MeV, 
i.e. near and below the particle thresholds, the deformation impact almost vanishes because of the 
compensation of the deformation contributions from /i = and = 1 GDR branches. The effect 
of triaxiality should still be checked. We plan to extend our exploration to Nd and Sm isotope 
chains where the triaxiality is absent. 
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